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1. INTRODUCTION 


The developments of new aircraft, space rockets and guided missiles 
require a vast amount of research, design and construction, the cost of 
which frequently runs Into millions of pounds and many months of time. 

Any processes that can reduce this expenditure of time and money is, 
of course, vigorously pursued. The wind tunnel testing method is a 
convenient method for this purpose. 

Most problems of experimental aerod}mamlcs are connected with the 
study of motion of a body in relation to a stationary fluid. However, 
we can reverse the problem and study the motion of a fluid In relation 
to a stationary body. 

Dynamical similarity Is the basis of wind tunnel testing and other 
forms of model testing. 

According to the theory of similarity, the experimental results from 
testing small model can be converted to the full scale object. Usually 
the tested model is mounted in the tunnel with the aid of different types 
of suspensions, supports and struts, etc. But these support systems will 
bring about a series of additional questions. Their influence on the flow 
pattern in the tunnel and around the tested model is considerable. In 
the general case, these effects are expressed in changes in the velocity 
and pressure distributions, which are noticed: 

a) As changes in the average velocity in the testing section which 
necessitates corrections in the velocity coefficient of the tunnel. 

b) As changes in the pressure gradient, which create a horizontal 
Archimedian force affecting the drag, then necessitating a correction in 
the pressure gradient. 

c) As different local influences affecting boundary- layer flow, 
vortex formation, local flow separation etc. Especially in a supersonic 
wind tunnel, any support members attached to the model ahead of its 
most rearward point will create a shock-wave disturbance and affect the 
flow aft of that point within a roughly conical region. 
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So it Is better to have a model which can be suspended in the wind 
tunnel without any strut. 

Aeronautical scientists have achieved this technique successfully 
by using an electromagnetic technique. 

A model which consists of a permanent magnet is placed in a wind 
tunnel around which there is a set of electromagnets. The permanent 
magnet is attracted by the electromagnets and thus can be suspended in 
the wind tunnel. The force that the electromagnets attract the permanent 
magnet is proportional to the magnetic moment per unit volume of the 
permanent magnet multiplied by the field gradient produced by electro- 
magnets . 

Because the magnetic moment of the permanent magnet is limited to 
the saturation value and large field gradients are limited by the size 
of aperture containing the tunnel and model, the size of the testing 
model is limited. 

To increase the field gradient would be prohibitively expensive 
whereas to increase the magnetic moment per unit volume of the suspended 
magnet appears to be feasible through the use of a superconducting magnet 
to replace the permanent magnet. In practice, a combination of super- 
conducting magnet and permanent magnet may be necessary for control 
purposes . 

In order to keep the superconducting magnet in operation, a 
horizontal dewar containing liquid helium is to be designed and built 
for the purpose of this project. 
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THE STRUCTURE OF THE HORIZONTAL DEWAR 


This Dewar is made up of the inner can 13, the outer can 15, the 
liquid helium filling tube 17 and the evaporated helium vapour vent 
tube 20. The superconducting magnet solenoid 12 and the superconducting 
magnet switch 16 are in the inner can. There is a vacuum space between 
the inner can and the outer can. This vacuum space is filled by the 
multi-layer super- insulation which consists of about 20-25 layers of 
Mylar and carbon- loaded paper. 

The inner can and the outer can are made up of thin-walled tubes of 
stainless steel because the combination of reasonable strength and high 
toughness offered by autenistiC stainless steels makes them suitable for 
use at low temperature. Another factor that must be considered is the 
weight of the Dewar. The thinner wall thickness of the stainless steel 
is suitable for reducing the weight of the Dewar. 

The spiral tubes are used both in the filling tube and the vent tube 
for increasing the length of the tube and thus minimizing the heat Inleak 
from heat conduction. Another function of the spiral tube is to maintain 
a horizontal liquid level when the Dewar is flying horizontally. 

In order to avoid vacuum collapse, stiffening bands of aluminium are 
used to wrap the outer can. 

There is a vacuum valve in the left-hand cap of the outer can. 
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1. Bush of -vacuum valve (Brass) 

2. Vacuum valve (Brass) 

3 . Outer Can Cap(Brass) 

4. Inner Can Support (Tuf non) 

5 . Nuts (Nylon) 

6 .Bolt (Nylon) 

7 .Bush (Brass) 


Fig.l. The 


S.Imer Can Cap (Brass) 

9 . Stiffening Band( Aluminium) 
10. Shield Foil (Aluminium) 

11. Former of Coil 

12. Magnet Coil 

13 . Inner Can(StainlessSteel) 
l4.Super-Insulation(Al-Foil) 


Structure of the Flying Bev/ar, 


15 . Outer Can (Stainless Steel) 

16. Super-conducting magnet Switch 
17 . Filling Tube(Stainless Steel) 
18 . Inner Can Cap (Brass) 

19 . Bush (Brass) 

20 . Vent Tube (Stainless Steel) 

21. Bush (Brass) 

22»0uter Can Cap (Brass) 
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3. CALCULATION OF THE HEAT INLEAK AND FINDING OUT THE VOLUME OF THE INNER 
CAN FOR OPERATION IN THE DESIRED TIMES. 

From the diagram of the Dewar, we note that the heat Inleak into 
helium Is mainly caused by 

1) heat conduction throught the filling tube, 

2) heat conduction through the vent tube 

3) heat conduction through the support rod and wheel. 

4) heat radiation through the superinsulation space, 

let us now consider these sources of heat Inleak Individually. 


3.1. Heat Conduction Through the Filling Tube 

When we deal with the heat conduction problems in material, we get 
the heat conduction differential equation (1) according to Fourier's 
Law and the Law of energy conservation; 


That is q 


„ 8T' 

^ S 8-t 


„ 3T 


V(KV T). 


(|^(k|^)+ ^ (kP) +r^ (k|^)) 

3x 3x 3y ,3y 3z 3z' 


..( 1 ) 

..( 2 ) 


Where is the heat generated in unit volume and in unit time. Joule/sec m 

p is the density Kg/m^ 

C is the specific heat Joule/Kgmole K 
P 

K is the thermal conductivity Joule. m/sec K. 


For steady state, the temperature distribution inside the solid will 
not change with the time. SO = 0 

If there is no heat source inside the solid, “ 0 
Thus the equation (1) reduces to V (K V T) = 0 


(3) 
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In the one-dimension problems the equation (3) reduces to 

3 3 


— (K — -) 
3 X 3 X 


(4) 


Consider a rod or a tube which is thermal insulated in the traverse 
dimension. The length of the rod is L, the two end temperatures are 
T. and T . A coordinate system is chosen as shown in Fig. 2. 

When K = constant, the equation (4) reduces to 


0 (5) 

The general solution of the differential equation (5) in the condition 


that X = 0, T = Tj and x = L, T = T^: 
T„ - T, 


X + T, 


According to Fourier's Law of conduction, Che heat flow rate 

Q 


-KA = — (T - T ) 
^3x L ^2'’- 


( 6 ) 


When K is dependent to temperature we can introduce a new variable E 

r 

o 


/■ T 

defined by the relationship E = ^(T) where T^ is a reference 

temperature. ° 

According to its definition we can get 

3 _ 3 , 3E 3t. 3^E 

3 X 3x^ ■ 3 X ^ 3t • 3x^ ~ 3 x' 


3 E 


so the equation (4) reduces to 


3.^E 

3x^ 


= 0 . 


(7) 


We can find Chat Che equation (7) is very similar to the equation (5). 
By similarity, we get the general solution of the equation (7) when 
X-0, E-E( J, X=L, E=E 


E = 


^(T,) ~ ^(Ti) 


KA 


X + E 


We also get: q = _ - E 


(Tj^) 

1 ... 


( 8 ) 
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Fig. 2. A rod with one dimension 
heat conduction 



Fig. 3. Heat conduction through the vent tube 
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Ue use the expression (6) to calculate heat flow rate In the case of one- 
dlmenslonal steady>state heat conduction where the thermal conductivity is 
considered to be constant. 

If the thermal conductivity of the material which makes up the 
conduction region varies significantly in the region, we use the expression 
( 8 ). 

At first let us outline our Dewar: 

The outer diameter of the spiral tube is ~ 0.64 cm. 

We have two turns of the spiral tubes. 

The outer diameter of the inner can is 2" ~ 5.08 cm. 

The outer diameter of the outer can is 2%" - 6.35 cm. 

The thickness of the coiled tube is 0.01" « 0.025 cm. 

The thickness of the inner can is 0.01" ■ 0.025 cm. 

The thickness of the outer can is O.Ol" •* 0.025 cm. 

Thus we can find out: 

The cross section area of the colled tube is 

A = rdt= rx( — ) X (2.54 x O.Ol) = 0.03067 cm^ 

The length of the coiled tube is 

L = ttd-2 = 7t X (2.54 x 2) x 2 = 31.9 cm = 30 cm 

To check the thermal conductivity of stainless steel, we get 

K = 3 X 10~^ watts/cmK at T = 4 K 
K = 150 X lO"^ watts /cmK at T => 300 K 

Assuming that there is a linear relation between the thermal conductivity 
of stainless steel and the temperatures. We find out the dependence of K 
to the temperature: 

K = Kg + A(T - Tp) = 3 X lo"^ + 0.5(T - 4) x 10~^ (9) 
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According to expression (8) 

(~300 

Q.L = = A(E3qqj^- E^j^) = A J ((K^ b(T-T^))dT 

« ( K^(300-4) + Y <300 - 4)^ ) A 
= ( 3 X 10“^ X 296 + ^ X 296^ x lO"^ ) A 
= 22.79A 

So Q = 22.79 -^ = 22.79 38.5 x lo"^ watts 

3.2. Heat conduction through the vent tube 

In this case the escaping vapour maintains good thermal contact 
with the wall of the tube and absorbs a great deal of the heat tha't would 
otherwise enter and cause evaporation of the liquid helium. 

For the control volume surrounded by the boundaries of cross section 
A, B and the tube wall, is the conducted heat which flows to the liquid 
helium from the wall. Q is the heat flux in the cross section A which 
flows into the control volume. We assume that the tube and the vapour 
in the same cross section have the same temperature. 

Obviously, the heat which is transferred into the control volume is 
equal the heat which flows out from the control volume. Thus: 

Q + mC T = W + m C T 
p o. o P 

So Q = m Cp (T - T^) + W^ (10) 

where m is the mass flow rate 

Cp is the specific heat of the vapour. 

According to Fourier's law of heat conduction, the heat flow in the 
section A is equal to 




(in the case of above coordinate system) 
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When K is dependent to the temperature, it is expressed by the expression (9). 

Then Q = A (K +a(T - T )) (11) 

o o. aX 

upon combining Equation (10) and (11) and integrating, it is found 
that: 

, W + m Cp A T 

L/A = ( ) (a A T + (K - W — ) In — ), . . (12) 

m Cp m Cp W^ 

where L is the length of the tube and A T the temperature difference 
between the warm and cold ends. 

In the present situation, Cp, the specific heat for helium at 1 atm 
is about 5.2 Joules/g.K. 

• A T = 300 - 4 = 296K 

-3 

K = 3 X 10 watts/cmK 

° -3 2 

a = 0.5 X 10 watts/cmK 

According to these figures, we can plot a graph of m against L/A for 

various W . It is shown in Fig 4. The calculated results are shown in 
o 

Table 1. 

At present, we can only estimate the mass flow rate m. We estimate 

_3 

the total inleak is about 200 x 10 watts. Thus the evaporated liquid 
helium by absorbing this amount of heat will be: 

200 X 10“^ 200 x 10“^ ,^-2 , 

= ss 10 g/sec 

L 20,42 

c 

where is the latent heat of liquid helium at 1 atm. 

From A = 0.05067 cm^ L = 30 cm we get L/A = 592. 

Now we can check the diagram m 'v, L/A for various W^. We find that 

W will be less than 10 Joule/sec. 
o 

So we can conclude that W^ is nearly equal to zero, W^ = 0. 

It means that the escaping vapour will absorb almost all the heat which 
would be transferred to the liquid helium by heat conduction along the 
vent tube. 
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Table 1. The calculated results of m against L/A for various W 

( 

according to Equation (12). 



hJ, ( t^cbtis ) . 

-to 

10 

10^ 

-s 

lO 


i6^ 

10^ 

10'' 

/ 


10-^ 

4o6 

353.6 

340,0 

325.6 

304,5 

249 

J/7.6 

20.7 

10^ 

4/.^ 

36.68 

3S3S 

34.0 

32. S6 

30.4 

24.9 

/A7& 


!0-^ 

4.33 

3.50 

5.63 

3.53 

5,40 

3.3-S 

3.£>4 

2.49 


I0‘' 

0Mb 

0.39 

0.530 

0.56 

0.355 

0.54 

0.324 

O.3o4 


/ 

0.046 

o.o40 

0.039 

o.m 

0.0567 

0.0553 

0054 

(f, 033.5 


3*3. Heat Conduction Through the Support Rod and Star-Wheel . 

The supporting system is composed of nylon rod and tufnol star-wheel. 

For the problems of heat conduction through the nylon rod we just 
simply use equation (6) but for the tufnol wheel the problem is- more 
complicated. 

At first, let us consider a long cylindrical element of uniform 
thermal conductivity, as shown in Fig 5. If the conduction process is 
in a steady state and if the conditions specified on each of the two 
surfaces are independent of 0 and z the process is one- dimensional. 

That is, the temperature distribution depends only upon y. If temperature 
of Tj^ and are specified at and y^, the relevant differential equation 
and boundary conditions are as follows: 



at Y = Yj^. T = T^ 
at Y = Y2> T = T^ 

The general solution is then T 

In determining the heat flow rate, we apply the Fourier equation in 
cylindrical form 





13 . 



Fig. 5 -The diagram of the cross section 
of a cylinder 



Fig. 6. A rod with varying cross section area. 
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= -KA(y) 


dT 

dy 


The area in this case is equal to 2 tty-I** (I* Is the length of the cylinder) 
Then we obtain 




2 ttKL 

in^ 


T - T 
2 

«2 - •'!> ■ - 


1 ^ 
2 ttKL Yi 


(13) 


We can express equation (13) as following form: 


T2 - 




2 wkL 


^ 2 

In = R cylinder 

Yi 


(14) 


T -T 
2 1 

Now we define R as the thermal resistance of the Cylinder. 

At = T 2 ~Tj^ is the te^erature difference in analogy with the electric 
voltage, is the heat flow rate in analogy with the electric current. 

Secondly, let us consider a rod of varying cross section area as 
shown in Fig 6. The boundary conditions are as follows, 
at X = 0, A = A^. T = T^ 
at X = L, A = A^. T = T^ 

The relation of the cross section area to X is a linear relation. So we 

A2“A3 

can obtain A^^^ = A^H — — X. The heat flow rate is obtained by applying 
Fourier Law of heat conduction. 


dT 


(^ is constant for steady state heat conduction) 


Separating variables and integrating, we obtain that; 




v-rod 


The total thermal resistance of six rods of varying cross section area 

as shown in Fig 7 is one-sixth of the R , 

° v-rod 
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so: 




R 


total rod 


(15) 


Finally we use equation (6) to get the thermal resistance of nylon bolt. 



nylon 


nylon 


K A 


nylon 


(16) 


Upon addition of equation (14), (15) and (16), the unknown temperatures 
(T^, Tj^) are eliminated. 


so £ = 


^ 3-^0 


total 


(17) 


Equation (17) is the expression for the heat flow through the supporting 
system in terms of the overall temperature difference. 

Let us p re- determine that: 

The length of the nylon bolt is 2.0cm 
The outer diameter of the bolt is 0.6 cm. 

The inner diameter of the bolt is 0.3 cm 

The outer diameter of the disc is 2 cm 

The inner diameter of the disc is 0.6 cm 

The thickness of the plastic plate is 0.2 cm 

For the rod of varying cross section area, its thickness is 2mm 


= 0.2cm X 0.2cm A 2 = 0.4cm x 0.2cm 
The length of each rod is about 20mm 

The thermal conductivity of plastic plate is considered to be 

_3 

constant, K = 3x10 watts/cmK 

The thermal conductivity of nylon bolt is considered to be constant 
K = 3xl0~^ watts/cmK. 


_ nylon 


2.0 


V» 3 K 10-3 j „ 32, 


= 3143 K/Watts 


cylinder 


In 




2ir 3 X 10 0.2 


-In — r = 319.36 K/watts 


2rr kl 
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v-rod 


. L 

'■ A^-A^ 6K 
In A 2 


So R 


total 


•^3 

- SR - 
n 


= In = 1925 K/watts 

(0. 4-0. 2)x0. 2x6x3x10 

3929 + 319.36 + 1925 - 5387.5 K/watts 


A T 300-4 o -irt-3 

So Q = n — “ - eTa ' T-c : “ 54.9 x 10 watts 

\otal 


3 . 4 . The Transfer of Heat by Radiation Through the Superlnsulatlon Space 

The major heat flow through an evacuated insulation is radiation. 

The fundamental equation for heat transfer by radiation from a 
black body (a perfect radiator) is the Stefan- Boltzmann Law. 

^ = 0 at"^ 

rate of emission of radiant energy 
surface area of the emitter 
surface temperature of the emitter 

-8 -2 -4 

Stefan- Boltzmann constant, 5.67.10 wm K 

Most bodies are not perfect radiators and emit radiation at a lower rate 
than a black body thus 

^ e A 

where e, the emissivity of the surface, is generally a function of T. 

If a body at a temperature equal to that of a black body emits a constant 
fraction of the black body emission at each wavelength the body is called 
a grey body and e is a constant for that body. 

If a black body area Aj^ at a temperature Tj^ is completely surrounded 
by another black body at a higher temperature then the net rate of 
radiant heat transfer to the inner body will be 

% = A^(t/- t^S 

If neither of the bodies is a black body and the one body does not 
necessarily surround the other then the above equation is modified as 
follows. 


where q = 
A = 
T = 
a = 
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ei 62 


62+ (^) (1-62)6^^ 


(T2^- 


If Aj^ is nearly equa? to and e " thermally floating 

reflecting surfaces with emissivity § the radiation heat transfer 
between parallel surfaces at temperatures and T 2 is given by 


ir 


n+1 


-2- (T T 

2-e ^ 2 h ■' 


An effective mean thermal conductivity K can be defined using the equation 

( 6 ). 


_S_ 

A 


T -T 
2 1 


(18) 


where t is the thickness of the insulation 
then 

T T ^ 

K = eg 2 "1 watts/mK 

n+1 2-e T 2 - 

It is known that K = 0.52 x 10 ^ watt/cmK for the multi-layer super- 
insulation which is composed by aluminium foil 12 and 12 Dexter paper. 

P P 

It is also known that, 

the length of the superconducting magnet is 16cm , 

We assume that the length of the inner can is (16+L)cm. 

Thus the surface area of the inner can is 


A - j d^ + (16+L) TTd 

“jx 5.08^ + (16+L)irx 5.08 = 40.5+255. 3+15. 96L 
“ 295.8 + 15.96L cm^ 

t = = A- - 35 - 5.08 Q 

2 2 

According to equation 18, 

T - T 

S = K 0.52 X 10-^ X (295.8 + 15. 96L) 


- (71.6 + 3.86L) X 10 watts 
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Now we summarize that 

-3 

the heat flow through the filling tube Is 38.5 x 10 watts, 

the heat flow through the vent tube is nearly zero, 

-3 

the heat flow through the support system is 54.9x10 watts 

-3 

the heat flow through the evacuated space is (71.6+3.86L)xlO watts. 
Thus the total heat leakage into the liquid helium from the surroundings is 

a “ ^ ili “ (38.5f54.9+71.6+3.86L)xl0"^ = (165+3. 86L)xlo"^ watts 



THE DETERMINA.TION OF THE VOLUME OF THE INNER CAN BY HEAT BALANCE 

4.1. The Volume of Liquid Helium Evaporated for One Milliwatt Heat Inleak 

Checking the table of thermophyslcal properties of Hellum-4 we got 
that at 4.22 K and one atmosphere, 

the enthalpy of the liquid helium = 9.711 (J/g) 
the enthalpy of the vapour of helium = 30.13 (j/g) 
the density of the liquid helium = 0.125 (g/cc) 
the density of the vapour of helium = 0.01689 (g/cc). 

Thus the latent heat of the helium at above states Is 
X = 30.13 - 9.711 = 20.42 (J/g) 

If the heat Inleak Is 1 milliwatt, the volume of liquid helium that could 
be evaporated Is: 

V » = ■ ■■■— = 0.3917 X 10 cm /sec 

X 20.42x0.125 

3 

that Is 1.410 cm /hr. 

For our design, we Intend to operate the superconducting magnet for half 
an hour. During that time, one milliwatt heat Inleak could make 0.705cc 
of liquid helium evaporate. 

4.2. Heat Balance 


Apart from the volume which Is occupied by the superconducting magnet 
the volume of the Inner can which will hold the liquid helium Is 

Vj^ = i L - ^ X 5.03^ X L = 19.87L (cm^) 

-3 

We have found out that the total heat Inleak Is (165+3. 86L)xl0 watts. 

For half an hour of operating time, the volume of the liquid helium which 

3 

would evaporate with that amount of heat Inleak Is (165+3. 86L)xO. 705 (cm ) 
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so we obtain that 

(165+3. 86L)xO. 705 = 19.87L 
solving this equation we obtain that 
L = 6.78 cm 

According to calculation, the length of Inner can will be 16+L =■ 22.78cm. 

In order to allow a margin, we shall consider the length of the Inner can 
to be 26cm. It means that we adopt L “ 10cm. 
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5. ESTIMATE THE TOTAL WEIGHT OF THE CRYOSTAT 
5.1. The Weight of Outer Can 

According to the length of Inner can, we detemine the length of outer 
can to be 39 cm, 

the thickness of the outer can t is 0.01" •= 0.025 cm, 
so the volume of the outer can is 

V = d.L . t = TT X 6.3 X 39 x(2. 54x0.01) = 19.6 cm^ 

>' 1 outer 

3 

the density of stainless steel is 10 g/cm 

so W „ , , = 19.6 X 10 = 196g 

stainless ° 

The volume of two caps of the outer can is 

= 2( 7 D^x t )+2(ir Dtl)= 2x(4- x 6.3^x 0. l)+2(rx6. 3x0. 1x0,5) 
cap 4 cap 4 

= 8,213 cm^ 

The weight of two caps is 

Wo = pV„ = 8.6 X 8,215 = 70.6(g) 

4 cap 

The total weight of outer can is 
W = Wj^ + W 2 = 70.6 + 196 = 266.6(g) 
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The total weight of the Inner can is 

W“W,,.+W = 105.4 + 48.5 = 153.9(g) 

shell caps ® 

5.3. The Weight of Stiffening Bands 

The thickness of each band is 3inin ( ^ = 3mm) 

The width of each band is 20mm ( b= 20mm) 

3 

These bands are made of aluminium of which the density is 2.67 g/cm 

3 

So the total weight of five bands is 2.67 g/cm 

Wc «= P.T- V., = P ,(5ifDlt) = 2.67x(5TTx 6.35x2x0.3)=. 159. 7g 
o • D AL S 1 

5.4. The Weight of Both Filling Tube and Vent Tube 

Either filling tube or vent tube is composed of colled part and 
straight part. 

i) Colled tubes of both filling tube and vent tube 
the diameter of the tube d = %" = 0.64 cm 

the length of each tube L = 30cm 

the thickness of the wall of the tube t = 0.01" = 0.025 cm 
so its weight is 

= 2x( dtL) P g Painless ~ 2x( rx x 0.01x2. 54x30 )xl0 = 30. 4g 

ii) The straight part of filling tube 

the diameter of it is d = 3/8" - 0.95 cm 
the length of it is L = 7.5cm 

the thickness of the wall is t = 0.01" ■ 0.025 cm 
so its weight is 

W- = dtl P = TT X 4 X 2.54x0.01x2.54x7.5x10 = 5.70g 
f .s 8 ® 

iii) The straight part of vent tube 

the diameter of it is d = %" = 0.64 cm 
the length of it is L = 6.0cm 

the thickness of the wall is t = 0.01" = 0.025 cm 
so its weight is 

W = x dtl P = Tf X % X 2.54 X O.Ol x 2.54 x 6.0 x 10 = 3.04 (g) 
v.s 
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5.5. Vacuum Valve System 
It Is made of brass. 

Its outer diameter dj^ = 1,2cm. The thickness of the wall t = 0.15cm 
Its length 1 = 1.5cm 
Its weight is 

W = w dtlP = w X 1.2 X 0.15 X 1.5 X 8.6 = 7.29 (g) 

For slide-block. Its outer diameter dj^' = 0.8cm. Its inner diameter 
d^' = 0.3cm. Its length X = 1cm. 

Its weight will be 

W = £ (d'2 _ P = 1*8-6 = 3.71 (g) 

so the weight of the vacuum system is 
W = 3.71 + 7.29 = 11(g) 


5.6. Bushes of the Cap of the Outer Can is about 5g. Bushes of the cap of 
inner can is also about 5g. 

Thus we obtain that the total weight of the whole cryostat is about 


W 


w + w + w 

outer inner stiffening 

can can bands 


W. u + 

tubes 


W 

v.s . 


+ W, 


bushes 


= 266.6 + 153.9 + 159.7 + (30.4 + 5.70 + 3.04) +11+10 


= 640.34 (g) 



25 . 


6. SETTING THE CENTRE OF GRAVITY OF THE CRYOSTAT NEAR THE CENTRE OF MAGNET 

According to the demand of operation we should set the centre of 
gravity of the Cryostat near the centre of gravity of the magnet as 
close as possible. 

Now the structure of the outer can and the inner can are fixed, the 
position of the inner can where it is inside the outer can is also fixed. 
So the only adjustment is the position of the magnet inside the inner can. 

Let us obtain the position of the centre of gravity of the Cryostat. 
We just simply use the Lever-rule. We assume that the centre of gravity 
of the Cryostat (B) is X cm away from the middle point (A) of the outer 
can and on the right of the latter. 

the distance of W valve is (202.5fX)mm 

the distance of W (inner can) is (25+X)mm 

the distance of W (coiled tubes )lsabout (125-X)mm 

the distance of W(straight tubes) is about (175-X)mm 

the distance of W(outer. can) is Xmm 

the distance of W (bushes of inner can) is (105-X)cm 

All forces and distances relative to the point B are shown in the Fig 8. 
According to the lever-rule, the total torque in clockwise direction to 
Che point B must be equal Co the total torque in anti-clockwise direction. 

So we get 

11 X (202. 5+X) + 153.9x(25+X) + (266.&+-159.7)xX 
= 30.4 X (125-X) + 8.74x(175-X) + 5x(195-X) + 5x(105-X) 

640. 34X = 754.5 
X “ 1 . 18mm 

So the centre of gravity of the Cryostat is l.lSmm away from the middle 
point of the outer can and to the right of the latter. 

Thus we can determine the position of the superconducting magnet 
inside the inner can. 

If we assume Che centre of gravity of magnet is in the middle point 
of the magnet that is 80mm away from Che left end of the magnet. 

As shown in Fig 9, the position of the magnet inside Che inner can 
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is 76.2 mm away from the left end of the Inner can or 23.8mm away from the 
right end of the Inner can. 
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THE CALCUIATION OF THE HELIUM GAS PRESSURE IN THE INNER CAN 

7.1. The Released Magnetic Energy 

After all the liquid helium in the inner can has evaporated, the 
temperature of the superconducting magnet will rise and exceed its trans- 
ition temperature. At that time it will be driven into normal state. 

Its resistance- less property will be lost. The electric current of the 
superconducting solenoid will decrease in a very short time. So does 
the magnetic field of the solenoid. During that time, the Induced emf 
due to the decrease of the flux inside the solenoid will drive the current. 
The stored magnetic energy in the space inside the solenoid will be 
released as a form of heat energy. It is known that the electric current 
inside the solenoid is 20A, and the inductance of the solenoid is 9 Henry. 
So the released heat energy is 

E=^Ll^=%x9x 20^ = 1800 Joules 

7.2 The Specific Heat of the Magnet 

As we have known, the winding wire of the magnet solenoid is 
composed of Niobium- Titanium filaments and the substrate of copper. 

Because copper has a much higher thermal conductivity than Niobium- 
Titanium, when the quench occurs, the heat generated can transfer to the 
liquid helium quickly and thus prevent damage to the operation of super- 
conductivity magnet. 

When we calculate the Increase of the temperatures of both super- 
conducting magnet and the helium gas, at first we must consider the 
specific heat of the superconducting magnet and of the helium gas. 

There are two contributions to the specific heat of a metal. Heat 
raises the temperature both of the crystal lattice and of the conduction 
electrons . 

According to solid-state physics, we can write the specific heat 
of a metal as 

“ ^latt + ^eL = T “ 234R(|)^ + rT 


(19) 
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where R Is Che Universal gas constant 
N is Avogadro's number 
Kgis Boltzmanr^'s constant 

Y is a measure of the density of the electron states at the Fermi 

surface 

9 is the Debye temperature. 

The specific heat of a metal in superconducting state is more 
complicated than in normal state. 

As we have known the superconducting state is one of the thermo- 
dynamic states of metal. When the temperatures of the metal decreases 
below the transition temperature in the absence of any applied magnetic 
field there will be a phase transition. The lowest energy state of the 
metal will become a superconducting state. This phase transition is known 
as a second-order phase transition. At the transition there is no latent 
heat, and there is a jump in the specific heat. It is true for type I 
superconductor that there is a sudden change in the specific heat as the 
metal goes from the superconducting into the normal state. (When the 
phase transition occurs in a constant applied magnetic field, the trans- 
ition is of the first order, and latent heat is absorbed). There are 
three thermodynamic states for type II superconductors, namely super- 
conducting state, mixed state and normal state. When the metal goes from 
the superconducting into the mixed state, there will be a large narrow 
peak in the curve of specific heat against T more like a " x -type" 
specific heat anomaly. When the metal goes from the mixed state into 
the normal state, there will be a sudden drop in the specific heat. (Both 
transitions occur in a constant applied magnetic field). 

For our calculation, we only consider the specific heat of the metal 
above the transition temperature. 

For simplicity of calculation, we assume that the winding wire is 
fully made of copper. 

1) At first, let us calculate the specific heat of copper at 4.2 K aiid 
at 20 K by using the expression (19) 
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Zt is known that 

Debye temperature of copper is 310 K 
the atomic weight of copper is 64 
Total weight of magnet is 0.85 kg 

So the number of moles of copper is 13.28 gmoles 

The specific heat of copper contributed by Lattice vibration at 4.2K is: 

234R(|)^ = 234 x 8/3143 x x lo"^ J/g mole K. 

The specific heat of copper contributed by the conducting elections 
is; (r' = 7.4 X 10 ^ j/mole K) 

C , = rT = 7.4xl0"^x 4.2 = 3.108 x lo"^ J/gmole K. 

V electron gas 

The specific heat of copper by lattice at 20K is 

C , = 234 R(i-)^ = 234 x 8.3143 x = 0-522 J/gmole K 

V iatt u jiu 

The specific heat of copper by electron at 20K is 

C , ^ = rT = 7.4 X 10“'^ X 20 = 0.0148 j/gmole K 

V electron 

So the heat capacity of copper at 4.2K is 

C = N(C , + C ,) = 13.28 X (4.838 + 3.108)x lo“^= 0.105 J/K 

V latt V el 

The heat capacity of copper at 20K is 

C = N(C , _ + C ,) = 13.28 X (0.522 + 0.0148) = 7.128 j/K 

V lat V el 

So we can see that even at 20K the heat capacity of copper is still 
rather small. For calculating the specific heat of a metal, the 
expression (19) is only applied in the low temperature region (T/6 < 0.1). 
So for higher temperature region we must use the diagram (10) to 
calculate the temperature rise of copper. 

2) Estimation of the temperature increase of copper by assuming that all 
the released heat energy is absorbed by the copper. 

We can see from the diagram (11) that either or E/T changes linearly 
with T/0 in the region T/0 = 0.1 '''0.3. 

Checking the diagram, we get: 
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3) Calculation of the final temperature of both copper and helium vapour. 
At first, we calculate the heat absorbed by the helium vapour and 
copper when the temperatures of the helium vapour and copper change 
from 4.2K to 31 K. 

Since the specific heat of copper In the temperature range of 31 K— ^ 

62 K changes linearly with T/0 we can get the average specific’ heat of 
copper In this range. The specific heat of helium vapour In this range 
Is about a constant. Thus we can obtain the average heat capacity of 
copper plus helium vapour. Then we calculate how many degrees the 
temperatures of copper and helium vapour can Increase due to the 
absorbed heat energy starting from 31K. At last we can obtain the 
final temperature. 

a) According to T = 31 K and P= 0.01689 g/cc, we get from the 
diagram of helium properties or from the table: 

= 108 j/g 

The absorbed heat energy by helium vapour when its temperature chances 
from 4.2K to 31K is:, 

Ae, ^ = (108-24.13) X 3.834 = 321.56 joules 

vdpour 

The absorbed heat energy by copper: 

when T/e =0.1 (T=31K), = 0.481 (Checking the table (II)) 

We consider the internal energy of copper at 4.2K as zero, so the 
absorbed heat energy A E is considered to be equal to E. 

A E = E = 0.481 T = 0.481x0.1x310 = 14.91 Joules/gmole 

A E = 14.91x13.28 = 198 joules 

copper ■' 

So the absorbed heat energy both by copper and helium vapour when 
the temperature of the system changes from 4.2 to 31K is 

^ = A E + A E = 198+321.56 = 519.5 joules 

total vapour copper 

b) The average specific heat of copper in the region (T/0=O. 1-^0. 167) 



35 


when T/0 = 0.1. | = 0.481. E 


14.91 Joule/gmole 


when T/0 = 0.167 | 
> T 

= — = 
^vaverage A T 


= 1.93. E = l,93x 0.167x310 = 99.92 Joule/gmole 

' ‘-O’ " 


The heat capacity of copper 


C = 4.09 X 13.28 = 54.32 joules/K 

copper •' 

The heat capacity of helium vapour Is 12 joules/K 

The total heat capacity Is 


^total 


c + c 

copper vapour 


12 + 54.315 = 66.32 joules/K 


c) The released heat energy Is 1800 Joules. 


When the temperature of the system goes from 4.2K to 31 K, both copper 
and helium vapour absorbed 519.5 joules of the heat energy. 


The heat energy to be absorbed Is A E = 
The Increase of the temperature will be 


A T = 


A E 


'total 


1280.5 ^ 
66.32 


1800-519.5 =■ 1280.5 joules 


So T = 31 + 19.3 => 50.3 K. 


d) Checking the calculated results. 

For T = 50.3 K. T/0 = = 0.162 

When T/0 = 0.167. | =1.93 

T/0 = 0.154 I = 1.59 

E 

By Interpolation, when T/0 = 0.162, — = 1.80 

So Ae^^ ^ = 1.80 X 0.162 X 310 x 13.28 = 1200 Joules 
copper 

Check the table of thermodynamic properties of helium, we get the 
Internal energy of the helium at T = 50. 3K and p = 0.01689 g/cc 
E = 167.6. 

^ ^vapour total “ (167.6 - 24.13)x 3.834 = 550 joules 
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A = 550 + 1200 = 1750 joules. 

There is only a very small error of the calculation. 

Thus we get the final temperatures of the magnet and the helium vapour 
to be about 50 K. 

7.3. The Pressure of the Helium Gas in the Inner Can 

We assume that the quench of the superconducting magnet is so fast 
that the helium gas which has absorbed the released heat energy and has 
been In the high pressure can't escape from the inner can during that very 
short time. 

So we consider that the density of the helium vapour In the whole 
process will not. change. 

According to T = 50. 3K and p = 0.01689 g/cc, checking "the 
thermophysical properties of helium - 4", we obtain; 

P = 18 atm. 


The Bursting Pressure of the Inner Can 


The equation for calculating the bursting pressure is 



where S is the tensile strength 

t Is the thickness of the tube 
D is the diameter of the tube 

It is known that, 

O O A 

the tensile strength of stainless steel is ISOxlO'^ Ib/ln ■■ 10555.5 kg/cm^ 
the thickness of the wall is 0.01" > 0.025 cm 
the diameter of the tube is 2" ~ 5.08 cm 


S = 150x10^ Ib/in^ = 10555.5 Kg/cm^ 


P = 2S.^ = 2^. 1 ^ 55.5 X 0.005 ^ atm. 

1.0336 


We adopt the safety factor as 5.5 
So the allowed working pressure is 
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working 


102 ^ 
5.5 


18.5 atm. 


So the pressure of the helium gas does not exceed the allowed working 
pressure. 

Thus we conclude that it is safe for the magnet to operate. 
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8. THE BASIC PRINCIPLES OF MAGNETIC SUSPENSION SYSTEM 


The principles of magnetic suspension system can be described 
by using the idea of the Interaction of the magnetic poles. It can be 
explained by the theory of magnetic field too. 

8.1. Description Based on the Magnetic Pole Theory 

According to classical magnetic theory, the force between two 
magnetic poles varies inversely as the square of the distance apart 
and is directly proportional to the 'strength' of the poles. This Is 
Coulomb ' s law . 

In order to suspend a model in a wind tunnel horizontally, a 
strong permanent magnet is mounted Inside the fuselage of the model. 

When the high speed air stream passes through the wind tunnel, there 
will be three directions of forces and three rotation moments acting 
on the model. (We adopt a system of orthogonal axes centred on the 
model, with one axis along the tunnel centre-line, a vertical axis, and 
a second horizontal axis transverse to the tunnel, say x,y,z axes res- 
pectively). They are x-direction force, y direction force, z direction 
force, moment about x axis, moment about y axis and moment about z axis. 
Hence the model will have transition movements in the x,y,z directions. 
The rotation about x axis called 'roll', the rotation about y axis called 
'yaw' and the rotation about z axis called 'pitch* also occur. In order 
to keep the model in equilibrium during the testing process, the magnetic 
suspension system should provide forces and moments to counteract these 
forces and moments. Furthermore those forces and moments acting on the 
model should be able to be measured. 

•Consider a pair of similar electro-magnets arranged around the wind 
tunnel as shown schematically in Fig. 11(a). Each magnet is separately 
energised. If the electro-magnets have equal energisation producing 
equal pole-strengths, the net force on a model positioned symmetrically 
below the electromagnets will act vertically through the centre 0 of 
the model. In this way the gravitational force and the Y direction 
force acting on the model by the stream can be counteracted. The net 
vertical force will be proportional to the sum of the currents in the 
colls of the two electromagnets. 
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Coil 1 Coil 2 



Fig. 11(a) Four force components acting on a model 


Coil 1 Coil 2 



Fig. 11(b) Resolution of forces into two components 
and a moment 


Coil 3 A Coil ^ 



Fig,.ii(c) Lateral magnet and drag magnet 
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F = K (N.I, + N,I,) 
y y 11 2 2 


where 


is a constant determined by experiment. 


If the electro-magnets are energised to different levels, the 
force on the model provided by the interaction between electromagnet 1 
and the model magnet will be different from the force of the interaction 
between electromagnet 2 and the model magnet. These forces can be 
resolved into vertical force F , horizontal force F and a moment about 

y * 

Z-axis see Fig. 11(b). The moment and the force F^ would be 
proportional to the difference between the strengths of the magnets 
(i.e. ampere-turns per metre of each magnet) 


M = K 
z mz 


(Njl, - DjI^) 


F' 

X 


fx 


(N2I2 - Njl^) 


where K and are constants determined by experiment, 
mz rx 

The force components in the y direction, as described before, would be 
proportional to the sum of the strengths of electromagnet 1 and 2. 

So we can control the force in y direction and the moment about 
z axis by changing the currents in electromagnet 1 and 2. 

(The X component force produced here is much smaller than the y 
direction force.) 

When the model moves away from its symmetrical position related to 
coll 1 and coll 2, a relatively large effect on the x direction force 
will be produced (see Appendix A). 

This can be used to resist the applied drag force if the drag force 
is rather small in the case of low speed wind tunnel. 

For a high speed wind tunnel, the electromagnet 7 is used to provide 
a force in the upstream direction (x direction) opposing the drag force. 
This force is proportional to the strength of the electromagnet 7. 


^x = Kf, 


N7I7 


where is constant d^t^ermlned by the experiment. 


41 . 


The lateral electromagnets are shown in Fig. 11(c). The diagram 
represents a horizontal section through the wind tunnel axis. A force 
component along the Z axis and a moment about y axis are provided by 
the electromagnets. Similarly the z direction force is proportional to 
the sum of the strengths of electromagnets 3 and 4 less the sum of the 
strengths of electromagnets 5 and 6. The moment about y axis is proportional 
to the sum of the strengths of electromagnets 5 and 4 less the sum of 
the strengths of the electromagnets 3 and 6. 

■ '‘fz »3'3 * ■'4'4 - Vi ■ W 


M 


my 


W515 + "414 - 


”3^3 




where and are constants determined by experiment. 

There are special arrangements for controlling the roll about 
X axis. Methods for roll control can be found in the reference material (1). 

8.2. The Description Based on the Theory of Magnetic Field 

1. Moment and force on a current loop 

A magnetic dipole experiences a torque in a uniform magnetic field 

> > > 

T = m X B 

where m is the magnetic dipole moment of a current loop, |m | = lA 

I is the current flowing in the current loop. 

A is the area surrounded by the current loop. 

1 is the Intensity of the uniform magnetic field. 

In a uniform field, the net force on a current loop is zero: 

^=^I(d£x^) = I_^dfxB =I(^d£)x^ =0; 

the constant ^ comes outside the Integral, and the net displacement 
^d£ around a closed loop vanishes. 
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In a non-uniform field, this is no longer the case. For example, 
suppose a circular wire of radius R, carrying a current I is suspended 
above a short solenoid in the "fringing" region (Fig. 12). Here ^ 
has a radial component, and there is a net downward force F on the 
loop (Fig 13), F = 2 tt IRB Cos 0. 

For an infinitesimal loop of dipole moment m in a field B 
the force is: 


^ = (^. V) I 

It is proved in reference (5”) 
It may be written in the form: 


F 


X 


> 9^ 

XU • ^ 

<?x 


F 

y 


> 3 ^ 

3y 



3B 
d z 


( 20 ) 


2. Force on a magnetic substance 

When placed in a magnetic field, a magnetic substance becomes magnet- 
ized. An element of volume dV of the magnetic substance acquires a 
magnetic moment. 

dm = ^ dV, 

where m is the magnetisation vector. According to equation (20), the 
force on this volume element of the magnetic substance is equal to 

dF = |-2 dv, dF = dv, dF = ^ dv (21) 

X 3x’y 3y’z 3z 


The integral of equation (21) over the whole volume of the magnetic 
substance will be the total force on that magnetic substance. 

3. Magnetic moments and forces on the model 

The actual magnetic moments or torques that can be applied to the 
model containing a magnetic substance or a permanent magnet are: 

^ ^ K V M X I (22) 

where K is a constant 

V is the volume of the magnetic substance 
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I 


is the average magnetization vector of the magnet, 
is the average applied magnetic field intensity. 


The magnetic forces are applied to the magnetized model by means of 
a set of controlled gradients in the magnetic fields. The general 
equation describing these magnetic forces is as follows: 

I «v- KV(^ . V ) 

This can be expanded as follows: 


F = KV 

X 

F = KV 

y 

F “ KV 
z 


(M ^ 
— X 3x 


+ M + M 
3x — z 


3B 3B 

(M + M + M 

— X 3y ~y 3y “z 

3B 3B 

(M ^ + M + M 

-TC 3 z —y 3 z — z 





(23) 


3b^ 

where 


Is the average value of 


3B 


X 

3y 


over model magnet volume. 


Others such as 



are similar as 



3B 

~z 

3 X 


3B 3B 3B 3B 

~z ~y ~z ~x 

3y 3y 3y 3z 


-Z2. and 
3 z 



By assuming the magnetization lies along the wind tunnel direction 
(that is ^ = o, ^ = o) we obtain 

y z 

3 B 


(24) 


F = KV Mx 5 -^ 

X — 3 X 

§_B 

F = KV Mx — - 

y - 3y 

3B^ 

F = KV Mx 
z — 3 z 

Thus, according to equation (24), clearly Independent control of 
the gradient of the magnetic field along and normal to the magnet- 
ization allows the force components to be controlled Independently, 
no matter what the direction of This is the approach that has 
been taken in the prototype system in fact. It is possible to 
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provide independent control of the three gradient components and 
at the same time maintain the magnetizing field in the vicinity of 
the model approximately constant. 

Equation (22) provides a means of estimating the magnetic torques 
generated through the Interaction of the magnetized body and the 
applied field. It is a ma-tter of little Importance whether ^ is 
due to a permanent magnet or is induced by a magnetizing field or 
is due to a equivalent superconducting magnet. The result is the 
same. Consider the first. From equation (22) in component form: 

= KV (My Bz - to By), 

T = KV (Mz Bx - Mx Bz) 

y - - - - 

= KV (to By - My Bx) 

Similarly assuming the magnetization is aligned with the x-axis, 

(I^= 0 , Mz = o) we also obtain: 

T = o 
X 

T = KV to Bz 

y - - 

T = KV to By 
z — —J 

In other words if the magnetization vector lies along the x or 
wind tunnel axis, a pitching moment or torque is developed by a 
vertical field and a yawing torque is developed by a lateral field. 
These two moments can be developed Independently as required. It 
means that these two moments can be controlled. 

The roll torque depends upon the difference in the curvature of 

B and B in the x direction. It is described in reference (3), 
y z 

Stephens and Goodyer Indicated the methods for generating an 
independent roll torque. The details are in the reference (1) and (2). 
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9. CALCULATION OF THE SUPERCONDUCTING MAGNET 


The structure of the superconducting magnet Is Just like a solenoid 
coll magnet, but the wire Is made of superconducting filamentary material 
N^- T^ In a cupronickel matrix. The length of the coll Is 2b » 120mm. 

Its outer diameter Is 2 ^2 °° 46mm. Its inner diameter Is 2 a^^= 18.9mm. 

The sketch of it is shown in Figjl4, 

NIOMAX CN A61/25 Is chosen for winding the coll. It Is an Intrinsically 
stable superconductor. It Is suitable for small colls without Inter-layer 
cooling. Its diameter Is 0.25 mm. If the external field is 6 tesla. Its ' 
critical current Is 27 amperes. The operating current Is about 21 amperes. 
The use of a twisted array of 61 nloblum-tltanium filaments" in a high 
resistivity matrix permits very rapid changes of current and magnetic 
field without quenching. 


The formula for calculating the central field is given as follows; 


(see ref. 4 )■ 


. (X + (OC 6^)^ 


H = J X a, 8 In 

O 1 


1 + (1 + 0 ^)^ 


where is the magnetic strength in central point (ampere/metre) 
is the inner radius of the coll, (m) 

0 is defined as 

^1 

a is defined as _2_ 

«1 

X Is the space factor. Is defined as 


^ ^ active section of the winding 

total section of the winding 
J is the current density in the conductor (ampere/m^) 


So we obtain: 


,-3 


= 9.45 X 10 m, a = 2.4338, 

’te.) . B i„ a+JaLtiiJ . 

1 + (1 + 


0 = 6.35 


1.378 
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The cross section of the wire is A = -r — 

. ird2 

The effective cross section area A* = 

2 

The number of turns of wires per m is n which equals 


^ 3 

ird^ irx (0.25x10"^)^ 


1.52 X 10 


According to definitions of J and 1 , 1 J = n I. 

where I is the operating current which equals 21.19 amperes. 


so 


B = Pq H = n I F(a, 6) 


= (4 irx 10“^) X 1.527 x 10^) x 21.19 x (9.45 x lo"^) x 1.378 
= 5.29 (Tesla) 


The volume of the wire winding space is 
V = 2 i: B(a^-l) a^^ 

= 2 TTX 6.35 X (2.4338^ - 1) x (9.45 x lo"^)^ 
= 1.65 X 10“^ m^ 


The total length of the wire in the coll is equal to 
J^= n.v. = 1.527 xlO^ x 1.65 x 10 ^ = 2521 m 


The number of turns of wires in the coll is equal to 
2 


N = 2aj^"(a-l) 6 n 

= 2 X (9. 45x10“^)^ X (2.4338-1) x 6.35 x 1.527 x 10^ = 24845 (turns) 


The Inductance of the coil is equal to 
, _ uoN^ A 

I 

Where A is the cross section area through which the flux of the magnetic 
field passes. 

2 

If we adopt A = ir 

. 4 TTX I0~^x 2.4845^x lo\ ttx (23xlO~^)^ 


120 X 10 


-3 


= 10.74 (henry) 
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If A is the average cross section area of the coll 

A = (a,+ a-)^ = 7^ X (23 + 9.45)\ lo"^ = 827 x lO'^ 

4 12 4 

L = 5.34 (henry) 

The results of calculation are summarized as follows; 

The central field Is 5.29 Tesla when operating current Is 21. 19(amperes) 
The total length of the wire is 2521(m) 

The number of turns of the wire is 24845 (turns) 

The inductance of the coil is 10.74(henry) 

The actual values of the superconducting magnet made by Oxford Instruments 
Ltd., are as follows: 

Central field is 5.3 Tesla 

Current for central field is 21.19 ampere 

Inductance of the coil is 9 henry. 
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10. THE EQUIVALENT AVERAGE MAGNETIC MOMENT PER UNIT VOLUME OF THE 
SUPERCONDUCTING MAGNET. 

10.1. Uniformly Magnetized Rod and Equivalent Air Field Solenoid 

It was Ampere's theory that the pronounced magnetic effects of an 
iron bar occur when large numbers of atomic-sized magnets associated 
with the iron atoms are oriented in the same direction so that their 
effects are additive. They may be regarded as tiny magnets or as 
miniature current loops. We can describe them by their magnetic moment, 
which can be expressed as lA'.A' is the area of a miniature current loop. 
I is the current flowing in the loop. The effect of the atomic magnets 
can be conveniently described by a quantity called the magnetization M, 
which is defined as the magnetic dipole moment per unit volume. 

A long bar magnet which consists of miniature current loops is 
shown in Fig. 15. The magnetic moment of each loop is lA'. 

Assuming that there are n loops in a single cross section of the 
rod. We have nA' = A, where A' = area of elemental loop. A = cross 
sectional area of rod. 

Further, let us assume that there are N such sets of loops 'In the 
length of the rod, then nN = N' 

where n = number of loops in a cross section or rod 
N = number of such sets of loops 
. N'= total number of loops in rod 
It follows that the magnetization M of the rod is given by 

ra _ N' I A' _ NI nA* NI 

“ ■ V LA T T “ IT “ 

where K'= equivalent sheet current density on outside surface of rod. 


L = Length of rod, m. 
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Fig. 15. Uniformly magnetized rod with elemental current loops. 
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Fig. 16. (a) Uniformly magnetized rod 
(b) Equivalent solenoid 
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From the end view of the rod, we can see that there are equal and 
oppositely directed currents wherever loops are adjacent, so that the 
currents have no net effect with the exception of the currents at the 
periphery of the rod. As a result there is the equivalent of a current 
sheet flowing around the rod. 

This type of a current sheet is just like the case of a solenoid 
with many turns of fine wire, as in Fig. 16. 

The actual sheet-current density K for the solenoid is 


K 


NI 

L 


(A m 


where N = number of turns in solenoid, 

I = current through each turn, A 
L = length of solenoid, m. 

If the solenoid of Fig. 16 is the same length and diameter as 
the rod of Fig. 16 , and K = K', the solenoid is the magnetic equivalent 
of the rod. 

Ve can obtain the magnetic field in the centre of the solenoid as 
follows (in the air) . 




At the centre of the rod, B = VJ<!= U M(Tesla) 

’ o o 

So we can also say that a uniformly magnetized bar magnet with 

magnetization M is the magnetic equivalent of the solenoid with a sheet 
NI 

current density K = — = M provided the rod is the same length and 
diameter as the solenoid. 


In this case, the total magnetic moment of the solenoid is 

NI 


M 


t , , . = NIA = . (LA) = K.V, 

solenoid L 


and total magnetic moment of the rod is 


M 


rod 


= M.V = K.V = M 


solenoid 
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So the total magnetic moment of the rod is the same as the 
solenoid which Is the magnetic equivalent of that rod. 

So If a solenoid Is presented. Its total magnetic moment divided by 
the whole volume will be the magnetization M of the equivalent bar magnet. 

10.2. The 'Average Magnetic Moment Per Unit Volume of the Superconducting 
Magnet . 

Our superconducting magnet Is Just like a multilayer solenoid. A 
representative current element 1 = JX dz dr Is shown In Fig. 17, where 
J Is the current density, X Is the filling factor. 

The area of the current loop which 1 flows through Is 
A = TT 

So the magnetic moment of the current loop Is 


dm=lA = Trr JX dzdr 


Integrating above equation over the entire cross section, we can obtain 
the total magnetic moment of the magnet. 


m = / 

o 


L iTni I* 2 

/ J A 7rr^ dr dZ = / / n I tt r dr dZ. (X J = n I) 


^1 '' ’^i 

So the average magnetic moment per unit volume of the superconducting 
magnet Is: 

li rm n *5 3 

/ n I irr^dr dZ n I(rm - 

- m o 

“ V 2 


u r m . L 


3 r^m 


It Is known that 


I = 21.19 ampere 

n = 1.527 X 10^ turns /m^ 

r = 9.45 X 10"\. r = 23xlo"^ m 
1 m 


So we obtain: 


_ p n I(rm^- r 

^ (Att xlO~^)x(1.527xlO^)x21.19x(23^-9.45^)xlO~ ^ 
3 X (23^ X 10“^) 


= 2.902 (Tesla) 
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So the magnetic equivalent of the superconducting magnet Is a permanent 
magnet with magnetization n^M = 2.902 Tesla. 

It Is not clear whether the magnetic moment of the superconducting 
magnet will vary with variations In the; gradient of the suspension 
fields. If It does, then position' control may be more difficult than with 
a permanent magnet. If the moment does not vary to a first order of 
approximation, then the control will be similar to that using a permanent 
magnet. Only an experimental test will reveal the answer. This will be 
carried out at a later date when the magnetic suspension system Is ready 
In Che Department of Aeronautics. 

If control proves to be difficult, the central hole will allow a 
high permeability ferro-magnet to be used for control, while the super- 
conducting magnet Is used for suspension. 
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11. THE EXPERIMENTAL RESULTS OF THE BOIL-OFF RATE OF LIQUID NITROGEN AND 
LIQUID HELIUM IN THE CRYOSTAT. 


The results are as follows: 



Boll-off gas flow 
(1 atm, room temp) 
(litre) 

Time required 
(second) 

Cryostat 

position 

liquid 

nitrogen 

0.5 

360 

vertical 

liquid 

nitrogen 

0.5 

340 

horizontal 

liquid 

helium 

5 

115 

vertical 

liquid 

helium 

5 


horizontal 


From above experimental results we can see that when the cryostat 
is in horizontal position the boil-off rate of liquid nitrogen or liquid 
helium is only slightly larger than that when it is in vertical position. 

When the cryostat is horizontal the temperature difference between 
the outer can cap and the inner can cap (that is between the ends of 
the spiral- tubes ) is equal to the room temperature minus liquid helium 
temperature or liquid nitrogen temperature. When the cryostat is vertical, 
the temperature difference between the ends of the spiral-tubes will be 
less than the former temperature difference because part of the stainless 
steel wall above the liquid level will become one part of the heat flow 
path and provide additional heat-flow resistance. Thus the heat Inleak 
is reduced. 

According to calculations from the experimental results, the heat 

-3 

Inleak with liquid helium in the cryostat is 144x10 watt. The super- 
conducting magnet can be maintained to operate at least 50 minutes. It 
has been proved in the experiment. 

So we can conclude that our cryostat is satisfactory for the require- 
ment of half an hour of operation time. 
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APPENDIX A 


It Is assumed that a model has a pole length S equal to the 
distance between the poles of the suspension electromagnets. It can 
only move horizontally In a plane beneath the electromagnets at a 
constant distance g, as shown In Fig. 18. 


It Is also assumed thatthe pole strength of the suspension magnet 
Is proportional to the current flowing In the coll of that magnet. Hence 
the forces between poles are given by 


IlP IjP V ^iP 

^1 - K 2 > ^2 “ ^ 2 ’ ^3 “ 2 ’ ^4 “ ^ 2 

“l ■ “3 ”4 

2 2 2 2 2 2 2 
where K Is a constant, m^^ = m2 = d + g , m^ =(s+d) .+ g . 

2 2 2 

m = (s-d) + g . P Is the magnetic strength of permanent magnet 
^ Inside the model. 

By taking vertical -components of f and f^, the vertical force on P^ which 
should oppose half weight of the model is equal to, 

II I, 

% y 

similarly the net vertical force on P2 Is equal to 
^2 ^1 

Kpg ( — 3 - — j) = %w 
m2 m^ 

Solving above two equations we get 


^2 = 


O'! 3 

s w(Mj^-" + ) 


2KP 


where M is non-dimensional form of m given by M = m/s, similarly 
G = g/s. 


and 


2 3 

3 w M, M, _ 

( + (_i )3 I • ) 

2KP ^ G ^ ^M3 ^ ^2 •' 


where I',= 12- Similarly one can define P, = — I- 

^ S w ^ 


S w 
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Both and the non- dimensional currents In the coll 1 and 
coll 2, are functions only of the geometry of the system. 

By solving the force components fj^ to f^ In a horizontal direction, 
a non-dlmenslonal axial force due to model displacement can be obtained 
In terms of non-dlmenslonal distance and currents. 


^ w Dl' DI-' (l-D)I ' (l+D)!.,' 

<— + + 3 ^) 

2 K, Mj-’ 


where D = — 
s 
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INSTITUTE OF CRYOGENICS 

THE UNIVERSITY 
SOUTHAMPTON 
S09 5NH 


Telex: 47661 

Telephone: (0703) 559122 
Extension: 


RESULTS OF EXPERIMENTS ON THE LIQUID HELIUM CRYOSTAT WITH A SUPER- 
CONDUCTING MAGNET PRIOR TO MAGNETIC LEVITATION 


We have carried out the following experiments: 

1. Cooldown of the magnet to the normal boiling point temperature 
of liquid nitrogen (77 K) 

2. Measuring the boil off of liquid nitrogen inside the cryostat, 
with the cryostat in vertical standing attitude and horizontal 
flying attitude respectively. 

3. Fill cryostat with liquid helium after precooling with liquid nitrogen. 

4. Measure the boil off liquid helium, with the cryostat in vertical 
standing attitude and horizontal flying attitude respectively. 

5. Energise the superconducting magnet and put into a persistent mode. 

Then allow magnet to quench as the liquid helium was totally evap- 
orated. 


1. The boil off of liquid nitrogen and liquid helium was: 



Boil off gas 
flow 

Litres at NTP 

Time required 
(second) 

Heat in 
leak mW 

Cryostat 

attitude 

Liquid 

nitrogen 

0.5 

360 

335 

vertical 

0.5 

340 

355 

horizontal 

liquid 

helium 

5 

120 

144 

vertical 

5 

108 

160 

horizontal 


It was noted that when the cryostat is in horizontal attitude the boil-off 
of liquid nitrogen or liquid helium is only slightly larger than that when 
it is in vertical attitude. 
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2. The boil off rates of liquid nitrogen when the vacuum space is 
pun 5 >ed for different times! 



Boil off gas flow 
(1 atm. room temp) 
(litre) 

Time required 
(second) 

Heat in 
leak mW 

pumping 
2 days 

0.5 

260 

465 

punp 
5 days 

0.5 

340 

355 


It was noted that the boil off of either liquid nitrogen or liquid helium 
is sensitive to the pumping time. After the cryostat is assembled, the boil 
off of either liquid nitrogen or liquid helium is lower, when the time to pump 
the vacuum space of the cryostat is longer. The reason is that the out- 
g^s sing process can be carried out more thoroughly via a long pumping time. 

In any case, we should maintain a vacuum in vacuum space and not expose it 
to atmosphere or we should have to pump it again for a long time. 

3. The time required for filling up the cryostat with liquid helium is 
80 minutes . 

The amount of liquid helium required in the transfer process is 4 litres. 

4. When we energise the magnet, the boil off is increased to 5 litres per 
minute of helium gas flow. This is twice as much as that in the normal 
operating process. 

It is suggested that we should energise the magnet before we stop the 
transfer of liquid helium. 

5. According to calculations from the experimental results, the heat inleak 
with liquid helium in the cryostat is about 150xl0“3 watt, equivalent to a 
boil-off rate of 215cc of liquid helium per hour. So we can conclude that the 
cryostat is satisfactory for the requirement of half an hour of operation time. 

6. The boil off of liquid helium with the cryostat at flying attitude is not 
affected by vibration of the cryostat. 

7. When we energise the magnet, the maximum current is 15 amperes. The 
maximum charging rate is 6 amperes per minute. 
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